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ABSTRACT 

Studies  have  been  published  which  extend 
and  apply  the  Kolomogorcv  theory  of  turbulence 
to  the  fluctuations  of  optical  refractive  index 
in  the  ocean  due  to  salinity  and  temperature 
nicrostructure . 

Mathematical  models  of  some  of  the  statis¬ 
tical  properties  of  these  fluctuations  are 
developed  to  provide  continuous  functions  obeying 
the  various  Fourier  transform  relationships 
provided  by  the  studies.  The  continuous  functions 
are  utilized  to  charecterize  the  phase  ana 
intensity  fluctuations  of  a  plane,  monochromatic 
optical  wave  pas.-ea  through  such  a  turbulent 
medium. 

The  knowledge  of  the  relationship  between 
fluctuations  of  the  optical  parameters  and  the 
refractive  index  micros tructure  provides  a  basis 
for  the  utilization  of  laser  beans  as  data  links 
or  in  studies  of  the  propagating  medium. 
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INTRODUCTION 


Considerable  effort  has  been  applied  to  descriptions  of  the 
temperature  and  salinity  microstructure  of  the  ocean.  These 
descriptions,  both  theoretically  and  empirically  derived,  are 
not  yet  definitive.  There  is,  however,  limited  agreement  between 
theory  and  measurement.  The  results  have  been  summarized  and 
discussed  by  Stanley. 1  The  descriptions  are  based  on  the  Fourier 
transform  related  structure  function  (with  distance  as  the  argu¬ 
ment)  and  spectrum  (with  wave  number  as  the  argument) .  These 
descriptions  are  composed  of  separate  functions,  applicable  in 
specific  ranges  of  the  arguments. 

For  a  similar  charecterization  on  a  measurement  plane  of 
the  optical  intensity  and  phase  fluctuations  of  a  plane,  mono¬ 
chromatic  wave  which  has  been  passed  a  given  distance  through  the 
turbulent  ocean,  the  descriptions  of  the  refractive  index  micro¬ 
structure  of  the  ocean  (based  on  the  temperature  and  salinity 
microstructure)  must  be  mathematically  smooth  to  eliminate  unreal 
effects  caused  by  "corners"  in  the  descriptions.  That  is,  a 
smooth  transition  function  must  be  found  between  the  descriptions 
in  the  separate  ranges. 

Expressions  valid  over  the  entire  range  of  interest  are 
developed  to  describe  the  refractive  index  fluctuations.  These 
expressions  are  shown  as  to  agree  with  present  theory  and  available 
data.  The  desired  descriptions  of  the  optical  fluctuations  are 
found.  As  some  of  the  statistical  descriptions  are  not  in  the 
common  literature,  the  analysis  techniques  for  obtaining  these 
descriptions  are  reviewed  in  appendix  A  for  both  scalar  quantities 
and  vector  fields. 


OPTICAL  INTENSITY  AND  PHASE  FLUCTUATIONS 


The  optical  intensity  and  phase  fluctuations  of  a  plane, 
monochromatic  wave  passed  through  a  turbulent  medium,  character¬ 
ized  by  refractive  index  microstructure,  will  be  described  on  a 
measurement  plane.  The  measurement  plane  is  perpendicular  to  the 
original  direction  of  cha  wave  and  situated  such  that  the  wave 
passes  a  distance  L  through  the  medium.  The  wave  number  of  the 
(light)  source  is  k,  the  wave  length  is  The  medium  is  described 
by  the  three-dimensional  spectrum  of  its  refractive  index  micro- 
structure,  Enn  (k0) . 


DESCRIPTIONS  OF  THE  INTENSITY  AND  PHASE  FLUCTUATIONS  BASED  ON 
SOLUTION  OF  THE  WAVE  EQUATION  BY  SPECTRAL  EXPANSION 


its* 


The  fluctuations  of  an  optical  wave  can  be  described  by 
amplitude  and  phase  fluctuations.  Ey  taking  the  ogarithm,  one 


^-Superscripts  refer  to  similarly  numbered  entries  in  the  Technical 
References  at  the  end  of  the  text. 
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separates  the  fluctuations  into  the  sum  of  the  logarithmic 
amplitude  or  intensity  and  phase  fluctuations,  which  can  then 
be  treated  separately  if  the  fluctuations  are  small.  Tatarski,2 
chapter  7,  applies  the  methods  of  "small"  and  "smooth"  pertur¬ 
bations  to  the  solution  of  the  wave  equation  by  spectral  expansion 
to  obtain  the  desired  statistical  descriptions  in  terms  of  the 
refractive  index  spectrum. 


The  two-dimensional  spectra  in  the  measurement  plane  of  the 
logarithmic  intensity  E£2^kp^  •  an<^  phase,  E^^p^'  are  given  as 


/  k  k2L 

E £2(kP’  =  **  L1 1 '  sin  ~ir 

P 

k  k2L 

1  +  — ry—  Sin  — 
k^L  k 


E.,(k  )  =  urt 

<p2  p 


Enn(kp) 


Enn (kp)  • 


(1) 


(2) 


The  corresponding  correlation  functions  are 


r 

£2<p)  "  2lTJ 

'  Jo<V,E£2<V Vkp 

(3) 

0 

r°° 

p2  (p)  =  2d J 

Jo(kpP)Ed2(kP)kpdkP 

(4) 

o 


The  corresponding  structure  functions  are 


°12  (p) 


V(p) 


Jo(kpP) 


E£2  (kp) kpdkp 


Jo  (kpP) 


E<?2(kp)kpdkp 


(5) 

(6) 


Appendix  A  reviews  these  statistical  descriptions. 

ASSUMPTIONS  REQUIRED  FOR  TK :  SOLUTION  OF  THE  WAVE  EQUATION  BY 
SPECTRAL  EXPANSION 

2 

As  stated,  Tatarski,  chapter  7,  applies  the  methods  of 
"small"  and  "smooth"  perturbations  to  the  solution  of  the  wave 
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equation  by  spectral  expansion.  These  methods  require  less 
restrictive  assumptions  than  those  required  if  the  principles 
of  geometric  optics  are  employed  (Tatar ski, 2  chapter  6). 

Tatarski  considers  a  medium  for  which  the  mean  refractive 
index  is  unity.  More  general  results  are  required  here,  but  the 
substitution  is  simple  and  does  not  affect  the  resulting  relation¬ 
ships.  The  derivation  of  the  relationships  involves  extensive 
work.  As  the  effort  here  merely  makes  use  of  the  results,  only 
a  brief  description  of  some  of  the  assumptions  will  be  given  to 
indicate  the  range  of  physical  validity  of  the  results.  The 
models  are  evaluated  beyond  the  range  of  physical  validity  only 
to  demonstrate  the  characteristics  of  the  functions.  The  assump¬ 
tions  are: 

e  The  wavelength  of  the  light  source  is  much  smaller 
than  the  smallest  Euclidian  dimension  of  the  inhomogeneities  in 
the  spatial  distribution  of  refractive  index,  tQ. 

X  <<  lQ.  (7) 

@  The  fluctuations  in  refractive  index,  n-.  ,  are  much 
smaller  than  the  mean  refractive  index,  nc,  where  tlie  refractive 
index,  n  =  n  +  n^. 

nl  ^  no‘  (8) 

©  The  wave  is  expressed  as  the  sum  of  a  perturbed 
component,  u^,  and  an  unperturbed  component,  uQ,  where  the  per¬ 
turbed  component  is  much  smaller  than  the  unperturbed  one.  More 
precisely,  the  assumption  is 

lul/u0l  ~  I n1/n0 |  <<  1.  (9) 


Let  log(u  )  =  ii0,  U]/u0  =  and  log(u)  =  Provided  that 
|ui/u0l  <<  1/  then  <!*]_  =  #  -  ij>0.  As  Inj/nl  <<  1,  |Vi|>|  <<  |^0I  is 
required.  As  |Vip0|2k  =  2tt/X,  it  is  required  that  Xl7tyj|  <<  2ir. 
Writing  'J'x  =  X1  +  iS]_,  it  can  tie  shown  that  |X]J  <  |Su,  so  the 
conditions  here  are  satisfied  if  the  phase  variations,  S^,  satisfy 


|  VSX  I  <<  27T/X.  (10) 


Under  these  conditions,  "the  angle  of  scattering  of  the 
waves  by  refractive  index  inhomogeneities  is  of  order  no  greater 
than  0Q  =  X/£q  and  is  thus  small.  Therefore,  the  value  of  ^  can 
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only  be  appreciably  affected  by  the  inhomogeneities  included  in 
a  cone  with  vertex  at  the  observation  point,  with  axis  directed  2 
towards  the  wave  source,  and  with  angular  aperture  9Q  =  A/f  <<  1." 

MATHEMATICAL  MODELING  OF  THE  REFRACTIVE  INDEX 
MICROSTRUCTURE  IN  THE  OCEAN 

In  order  to  permit  computer  aided  calculation  based  on 
equations  (1)  through  (6) ,  and  to  prevent  excessive  errors  due  to 
artificial  sharp  corners  in  Enn,  a  smooth,  continuous  function 
is  sought  as  the  mathematical  model.  The  difficulty  is  that  the 
available  theoretical  and  empirical  information  concerns  itself 
with  the  inertial,  viscous-convective,  and  decay  regions  of  turbu¬ 
lence  and  turbulent  mixing  separately.  Thus,  there  are  descrip¬ 
tions  of  the  refractive  index  variations  in  each  of  these  regions 
but  there  is  no  definitive  information  as  to  the  manner  in  which 
the  phenomena  behave  in  the  transitions  between  regions.  The 
criterion  of  selection  for  the  mathematical  model  of  the  spectrum 
is  that  it  and  its  Fourier  transform,  the  correlation  function, 
agrees  with  the  available  information  in  the  known  regions.  Once 
this  is  assured,  the  model  can  be  applied  to  describing  the 
optical  fluctuations. 

2  3  1 

On  the  bases  of  work  by  Tatarski  and  Laster,  Stanley  shows 
that  the  descriptions  of  refractive  index  fluctuations  (correla¬ 
tion,  structure,  or  spectral  functions)  are  a  linear  combination 
of  those  of  temperature  and  salinity  fluctuations.  These  func¬ 
tions  and  the  coefficients  of  the  linear  combination  are  scaled 
by  temperature,  salinity,  and  mean  refractive  index,  as  given  by 
Laster.  They  will  be  discussed  here  in  an  unsealed,  separated 
form. 

SPECTRAL  INFORMATION 

Table  1  gives  the  form  of  the  spectrum  in  the  three  regions. 


TABLE  1 
SPECTRA  FORMS 


Region 

One-Dimensional 

Three-Dimensional 

Inertial 

K«-VV5/3 

‘r  3  1 

Viscous- 

Convective 

-q*vl/2e_1/2k"*l 

fq^-V/V^V3 

i^vy/Vi/V3. 

Decay 

expi-D3/2k2q*/e1/2) 

exp(-D3/2k2q*/e1/2) 

a 
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In  the  inertial  region,  the  one-dimensional  spectrum  is 
taken  from  Gibson  and  Schwartz^  except  for  the  factor  1/2  which 
is  employed  so  that  che  one-  dimensional  spectrum  is  defined  over 
negative  and  positive  wave  numbers,  instead  of  just  positive. 

The  three-dimensional  spectrum  is  defined  based  on  equation  (11) 
from  Tatarski . 2 


1  3Eai(k) 
Eaa(k)  "  "  ~ 3k - 


(ID 


This  is  in  agreement  with  Stanley,  as  his  modified  three- 
dimension  spectrum,  has  a  factor  B2  =  |bt  which  corresponds  to 
the  results  of  equation  (12)  from  Hinz.e.5 


Ea(k)  =  4irk  £a(j  (k) . 


(12) 


£ 

Bp  is  a  universal  constant  given  by  Gibson  and  Schwartz'  as  0.31. 
The  factor  5/3  is  apparent  from  equation  (11) . 


A  similar  transformation  is  employed  in  the  viscous-convective 
region.  The  factor  q*,  given  as  approximately  2.0  by  Stanley-1-,  is 
given  by  Grant,  et  al,6  as  nominally  twice  that  (3.7±1.5).  The 
difference  is  apparently  due  to  the  difference  in  definition  of 
the  spectra. 


The  spectral  representations  in  the  decay  region  only  approx¬ 
imately  follow  equation  (11) .  The  form  of  the  three-dimensional 
term  is  found  in  Stanley^-  and  is  the  only  available  representation. 


7 

Batchelor  predicts  that  the  inertial  subrange  asymptote  and 
the  viscous  subrange  asymptote  (the  relations  of  table  1)  intersect 
at 


(13) 


The  intersection  calculated  for  one-dimensional  spectra  gives 


(14) 


The  intersection  calculated  for  three-dimensional  spectra  gives 
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2.15k  . 
S 


U5) 


k 


(3) 


1/4 


This  difference  has  occurred  in  the  literature.  Grant,  et  alf^ 
note  a  factor  of  (2.4+0  8)  x  10-2]Cs  working  with  a  one-dimensional 
spectrum,  while  Nye  and  Brodkey®  noted  a  value  5.0  x  10  2ks 
working  with  the  three-dimensional  spectrum. 

For  the  mathematical  model,  it  is  convenient  to  normalize 
the  one-dimensional  spectra  by  tne  factor 


q,vV2£-l/2/4lt  _ 


(16) 


making  the  one-dimensional  spectrum  in  the  viscous-convective 
region  simply 


2«k 


-1 


(17) 


The  result  for  the  inertial  region  is 


2-niB^q*)  k 


2/3t,-5/3 
s  " 


(18) 


STRUCTURE  FUNCTION  INFORMATION 

The  structure  function  information  is  taken  from  Batchelor7 
and  appears  in  table  2.  Of  course,  the  structure  function  rises 
to  a  constant  beyond  the  inertial  region.  In  the  inertial  region 
relationship,  a  proportionality  constant  is  needed  which  Batchelor 
assumes  is  unity. 


TABLE  2 

STRUCTURE  FUNCTION  FORMS 
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SELECTION  OF  A  MODEL 

Representative  values  of  the  known  parameters  are  as  follows: 

-3  ]_ 

£  =10  kinetic  energy  dissipation  rate 

—2  3 

v  =10  kinematic  viscosity 

a 

=  0.31  universal  constant. 

This  qives  ks  =  5.623  from  equation  (13).  Selecting  q*  = 
3.815,  k'~'  =  0.023ks  from  equation  (14)  and  k^)  =  o.049ks  from 
equation  (15) .  Let 


k3  =  k{3)  =  0.28; 


therefore 


k(1)  =  k3/2.15. 


3  —3 

where  for  representative  values  for  temperature  D  =  Dy  =  10 

(^5T  =  J-8)  aiid  for  salinity,  D  =  Ds  =  10-5  1x53  =  ISO). 

Equation  (18),  for  the  inertial  region,  may  now  he  written 


2-(3/5)k3/3k  5/J, 


and  the  viscous-convective  decay  regions  may  be  combined  as 
equation  (20) . 


2~k~1  exp  [-2 (k/k.) 2] 


Equation  (18)  for  the  inertial  regicn  and  equation  (20)  for 
the  viscous-convective  and  decay  regions,  when  considered  in  their 
respective  ranges,  provide  the  same  information  as  table  1;  that 
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is,  a  one-dimensional  spectrum  varying  as  k  ^  until  113/2.15, 
then  as  k~^  until  approximately  k.5,  then  an  exponential  decay. 

The  viscous-convective  and  decay  regions  are  combined  as  k3  is 
small  compared  to  k$. 

Permitting  the  spectrum  to  vary  as  all  the  'ay  back  to 

k  =  0  is  equivalent  to  assuming  an  infinite  energy  source,  wnich 
is  not  physically  realisable.  Instead,  equation  (19)  will  be 
said  to  hold  only  for  k  larger  than  kj,  where  kj  is  much  smaller 
than  >13.  This  provides  for  what  is  elsentiallv^a  scaling  of  the- 
problem.  The  selection  of  a  value  for  ki  sets  the  variance 
(the  value  of  the  correlation  function  at  zero  argument)  and  the 
upper  constant  of  the  structure  function  (the  value  at  infinite 
argument  which  is  twice  the  variance) .  For  k  less  than  ki ,  the 
spectrum  is  a  constant. 

The  forms  of  the  spectral  and  structure  functions  are  sum¬ 
marized  in  table  3.  The  structure  functions  are  normalized  by 
e juation  (16)  as  were  the  spectra. 

TABLE  3 

SUMMARY  OF  STRUCTURE  AND  SPECTRAL  FUMTION  FORMS 


lT.sxfS  !.ur£bcr 
?4Rce 

Spectres 

Structure  Function 

Distance 

Hange 

*.  », 

■■■  ■  -■  1 

Constant 

Constant 

r  >  1/k, 

ki  *  k  *  k3 

2- {3/S  ).-.%/3k~5/S 

2r(12.3) (3/S) (k,r)2/3 

l/kj  '  r  *  I  /': , 

h  *  k  '  k5 

2-*'1 

2-  lnOc*rV-.) 

1  rA  <  r  <  1/k- 

I 

! 

P 

f 

■*  “5 

*  o*p[-2fk/ksj‘ej 

2-11/6) fk^r)2 

r  *  1/k- 

1 

j 

Equation  (21)  is  a  smooth  function  whic‘  satisfies  table  3 
for  the  spectrum. 


(3/5) k^73 


v2 

“i 


1 5/6 


2  2 
5„  " 


a/2 


exp  -2 (k/k5) 


(21) 


Tatarskx  gives  the  following  Fourier  transform  pair  and 
related  information. 
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B(X)  =  /  COS(C5T)t:f-j)d» 


B  (t)  = 


- 1 H^L -  .  %  T) 

2a  XF (a  +  1/2)  f  a  ° 


B(0)  =  rliTT 


3(®>  -  0 


W(u)  =  2-  - 


2 \a+l/2 


= 


4^/T  L 

:a  1/2)  I* 


(~oT) 

(a) - ; —  S  (cs  t ) 

2°  1  a  °  . 


D<“i  -  23 to)  =  ,l°y t2) 


0(T)  = 


teVT  _i_ 

:(a  *r*  1/2)  2a 


r  (a  -  1>  (u:oT)^a 


;  ^ct  <  1, 

Eat  <  I 


where  B  is  a  correlation  function,  D  is  the  associated  structure 
function,  and  W  is  the  associated  spectrum.  As  the  first  fcerr. 
in  the  bracket  in  equation  (21)  ceccr.es  negligibly  srall  before 
the  exponential  multiplier  departs  significantly  from  unity,  the 
exponential  can  be  neglected,  bribing  this  term  as 


2  5/6 


and  considering  equation  (26)  with  „0  =  kj  .  a  =  1/3.  u  =  J< ,  we 
can  see  that  the  corresponding  structure  function  can  be  derived 
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■>  v  .! ..  k  '.r1-  I-’;  kJ&S)*nh" i  i- .fidii".  i*. 


from  equation  (27) .  At  small  values  of  the  argument,  equation 
i':9)  is  valid,  and  tb”  structure  function  is 


\3  . 


4  dl 

r(5/6) 


2/3-1 

r(2/3)  (k.r) z/-3 


=  2m (48. 25)  (3/5)  G^r)^7;  r  *  l/k3  . 


By  performing  an  approximate  integration  for  the  second 
tern  in  equation  (21),  including  the  multiplicative  exponential, 
one  gets 


sin  k-r 

8-  1*  In  kl  -  In  k, - —  -f  Ci(k,r) 

3  3  ,-3r  j 


-  cic:<:r)] 

•*  a 


-  8- 


j(-l)n:I{k3r)2n  (-l)nr2n(k2n  -  k'2n)l 


(2n  *  1) i 


2n(2n) ! 


wnere 


2- (k£r) 


r  v  l/k_ 


[o.577  -t-  In  (k*r)]  ;  l/k«-  <  r  <  l/k3 

I  e-[l  +  In  (k^/kj)]  ;  r>lA, 


r 

_  _  *  cos  t 

\4>#  “  *  ""“"T" -  Gu 


y<-»vn 

A2rt2n,! 


zr  r.e  cosine  integral  function  and  v  =  0.577  is  Euler's  constant 

The  factor  kj  is  the  wave  number  in  the  neighborhood  of  5:5 
at  which  the  exponential  decay  can  be  considered  to  begin  to 
strongly  affect  the  function  for  the  purposes  of  the  approximate 
integration.  If  is  selected  so  that  in  the  range  r  <  i/k$r 
the  result  in  equation  (32)  matches  the  desired  results  in  table 
3 ,  then 
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(33) 


(kji/kj-)  =  ~\Jl/5  =  0.408. 


This  is  a  reasonable  selection  for  k 5#  as  the  exponential  multi¬ 
plier  in  the  spectrum  evaluated  at  k£  is  0.72,  and  it  decreases 
rapidly  thereafter. 

The  desired  result  for  r  <  1/k-  no  longer  exists  when  the 
results  in  equations  (31)  and  (32)  are  added  for  r  <  I/k3.  In 
order  to  eliminate  these  undesired  effects  of  the  r  '  J  dependence, 
the  spectrum  model  is  changed  to 


(3/5)k2/3 


(3/5)\?/3 


(k2  *  k2)5/6  (k2  k2) 


ik2  +  k2)1/2 


,xp  [-2 


2I 


exp  I -2 (k/kc)  I 
t  3 


(34) 


In  finding  the  structure  functi.cn  associated  with  the  first 
two  terms  of  equation  (34) ,  the  exponential  can  be  neglected. 

The  result  is 


(Dr 


51  F (2/3) 


r (1/3) 


(k3/kx) 


2/3 


22/3  tor  r)l/3 

£  I  ***3/  ■*!  *  i'q£  / 


KV3,iclr' 


(k3r)1/3Kl/3(k3r) 


/  c 


2- (21.55)  (3/1'  (k,r) 


31.74 


o  ^  ,2/3 


r  <  1/i 


2/3  -  31.74  ;  1/k,  <  r  <  1/fc 


;  r  >  1/k, 


(35) 
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Combining  this  with  the  result  obtained  fron  the  last  tern 
of  equation  (34)  as  given  by  equation  (32) ,  we  get. 


f  31-74  [  (k3A1)2/3  -  l]  ; 


D„l<r> 


= 


l/kx 
2/3] 


l s[-  33  -  In  (x5A3)  +  (21.55)  (3/51  (k3r)"/3 


2-[o.IG3  ln(k*r4/4)]  ; 


1/k,  <  r  <  l/ke 


1/k-  <  r  <  1/k 


2- (1/6)  (kj-r)‘ 


Ccr paring  the  results  i 
on  table  3  we  see: 


r  <  i/k_ 


(36) 


n  equation  (36)  to  the  information 


•  The  constant  result  for  r  >  l/k1 . 

m  The  desired  two-thirds  dependence  for  1/k-  <  r  <  i/ki . 
This  ter-  is  multiplied  by  approximately  1.75;  but  as'stated  above 
table  2,  a  proportionality  constant  of  unity  was  assured,  so  that 
1.75  is  not  unreasonable. 


©  The  cesxrec  xogar: 
-ealicible  additive  constant 


~ic  tern  for  l/ke  <  r  <  1/k 


#  Tr.e  e-esxrea  scuare  lav  cepenoence  tor  r 


1/k, 


T«us.  tne  sexeccea  acei  rcr  tne  cne-axcensxcnax  si 
in  equation  (34)  Beets  the  specified  criterion  that  its  structure 
function  agree  with  existing  data.  Kote  also  that  the  intercepts 
of  the  asymptotes  of  the  co-del  are  still  as  stated  in  equations 
(14)  and  (15).  As  shown  in  table  1,  the  three-dimensional  model 
does  not  exactly  correspond  to  the  cne-dimensicnal  co-del.  In  the 
cede I  given  below,  an  error  term  exists  which  can  be  seen  to  be  a 
multiplicative  factor  on  each  term  consisting  of  unity  plus  a 
constant  on  the  order  of  unity  times  the  factor  %k/k^1~.  The 
factor  has  a  negligibly  small  effect  except  in 
The  exponential  for  the  decay  region 
the  references  to  cut  off  what  would 
biitlOB  fcn©  h~l  *  In  l 


is  heuristicaliy  chosen 
be  an  infinite  power  cor 
his  lioht,  consider ir.g  t 


tne  error  term  trees 
exponantial  itself  j 
is  neolioible. 


exponential 


ii  i 


tne  cecay  region, 


s  on  as 
he  aenroj 


ipidiy  as  ii 
alien  erres 
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grals  given  in  equations  {5)  and  (6).  A  computer  program  was 
prepared  by  Dr.  F.  Thsilheimer  and  Mr.  R.  Kybraniec  of  the 
Computation  and  Mathematics  Department,  Bethesaa,  to  perform 
these  integrations  by  means  of  the  eight-point  Gauss  quadrature 
formula.  A  sample  program  listing  and  additional  information  are 
given  in  appendix  B. 

Each  integral  is  calculated  with  unequal  intervals  starting 
at  0  and  endino  at  the  first  interval  past  1000.  The  intervals 
are  of  increasing  length  because  the  greatest  degree  of  fluctu¬ 
ation  occurs  for  the  small  limits  of  integration. 

Figures  5  and  6  display  the  structure  functions  of  the  phase 
variations  for  temperature  and  salinity,  respectively.  Figures 
7  and  8  display  the  structure  functions  of  the  logarithmic 
amplitude  variations  for  temperature  and  salinity,  respectively. 
Figure  9  displays  the  variance  of  the  phase  variations  for  both 
temperature  and  salinity.  Figure  10  displays  a  similar  result 
for  the  logarithmic  amplitude  variations. 

DISCUSSION 

The  mathematical  model,  equations  (37)  through  (39) ,  is 
shown  to  be  consistent  with  the  structure  and  spectral  function 
information  based  on  the  referenced  analytical  and  experimental 
results.  The  functional  form  of  the  node I  is  such  that  computer 
aided  computation  of  the  optical  propagation  properties  is 
possible.  The  smoothness  of  these  computer  results  partially 
demonstrates  the  suitability  of  the  model  and  the  computational 
technique. 

The  primary  results  are  given  in  figures  5  through  10. 

These  are  the  results  of  the  integrals  given  by  equations  {5}  and 
(6).  These  integrals  are  two-dimensional  Fourier  transforms. 

For  purposes  of  discussion,  the  term  integrand  will  refer  to  the 
two-dimensional  spectra  given  in  equations  (1)  and  (2) ,  each 
multiplied  by  the  variable  of  integration  k  .  The  functions 


will  be  referred  to  as  the 


wxnaow. 


Tne  results 


1 


<  0.1 


(40) 


(41) 


10 


will  be  utilized.  Thus  the  window  for  the  logarithmic  amplitude 
variations  is 


kpL 


V2r 

V* 


<  0.1 


>  10 


(42) 


while,  for  phase  variations,  the  window  is 


2 


< 


0.1 


k“L 


>  10. 


(43) 


As  the  three-dimensional  spectrum  has  an  insignificant 
contribution  to  the  integrals  for  k  <  ,  the  window  does  not 

affect  the  integrals  if 


10 


(44) 


L  > 


10.< 


1.32  x  1012. 


(45) 


The  term  "L  is  large”  or  "large  L"  shall  refer  to  L  >  1.32  x  10“  . 
When  L  is  large,  as  the  window  does  not  have  effect,  L  linearly 
scales  the  results,  but  has  no  affect  on  the  shape  of  the  struc¬ 
ture  function.  The  factor 


1  "  Jo(kpp)  l46) 

* 

shall  be  referred  to  as  the  transformation  function,  and  the 
result 
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0; 


1 


will  be  employed. 


Jo<V>  =  < 


1; 


kpp  <  0.1 

k  p  >  10 
P 


(47) 


The  results  will  be  discussed  first  in  terms  of  the  functional 
relationship  to  the  separation  distance  in  the  measurement  plane, 
p,  then  as  a  function  of.  L,  the  distance  traversed  through  the 
medium. 


When  L  is  large,  the  integrand  contributes  no  significant 
"energy"  (or,  to  paraphrase  Batchelor, ?  no  temperature  or  salinity 
stuff)  for  k  <  kj,  thus,  the  structure  function  should  stop 
increasing  for  p  <  1/k^  =  103,  as  the  results  show  (see  the 
l  =  1013  curves  for  this  portion  of. the  discussion).  As  p  is 
decreased,  the  transformation  function  inhibits  the  contribution 
of  the  integrand  to  wave  numbers  on  the  order  of  1/p  and  above. 

At  kc  the  integrand  decays  on  its  own.  As  p  decreases  from  1/k-, 
to  l/k3  (the  inertial  range) ,  the  results  show  a  p3'2  dependence. 


the  decay  range. 


a  p3  dependence  is  seen. 


|^arski' 


For  p  <  l/k5, 

predicts  the  pz  dependence  in  the  decay  ranee  but  finds  p' 
in  the  inertial  .range.  His  result  is  based  on  approximate  ana¬ 
lytical  integration  techniques  that  are  not  experimentally  con¬ 
firmed.  No  further  discussion  is  possible,  based  on  existing 
information,  except  to  otate  that  his  results  can  be  duplicated 
by  the  computer  techniques  utilized  here,  and  by  analytical  tech¬ 
niques,  assuming,  as  he  did,  the  inertial  range  spectrum  asymptote 
(k q-11/3)  exists  for  all  k  ,  rather  than  utilizing  the  given 
model. 


In  the  viscous-convective  range,  l/k5  <  p  <  l/k3,  the  results 
show  a  rise,  that  is,  the  curve  is  above  the  inertial  and  decay 
asymptotes.  As  this  exists  in  a  range  of  p  ir  which  physical 
measurements  can  be  made,  the  phenomenon  is  significant.  As 
Tatarski2  worked  with  a  medium  with  no  viscous-convective  range, 
he  had  no  such  results.  Similarly,  as  his  spectrum  rises  to  infin¬ 
ity  as  k  approaches  zero,  his  results  show  no  flattening  in  the 
structure  function.  The  f latte'  r.j  represents  the  finite  "energy" 
driving  the  turbulence.  If  a  smaller  value  of  k^  was  selected, 
say  at  the  scale  size  of  the  ocean,  it  is  apparent  that  the  variance 
would  be  larger.  The  curves  are  easily  extended  for  any  k^  less 
than  the  selected  1000  reciprocal  centimeters,  and  it  is  clear 
that  the  selected  value  of  k^  was  small  enough  not  to  affect  the 
results,  that  is,  to  permit  extrapolation. 

The  L  dependence  of  the  shape  of  the  phase  structure  function 
is  slight.  The  shape  of  the  curves  is  not  a  function  of  L,  while 
the  amplitude  is  linearly  dependent  upon  L.  Figure  9  shows  that 
at  L  of  the  order  of  1013-,  the  coefficient  of  the  linear  dependence 
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of  the  variance  on  L  ai\  *s.  This  is  noted  in  figures  5  and  6 
also.  At  large  p,  a  sep  .ation  change  between  curves  occurs 
at  large  L.  As  p  is  decreased,  the  separation  change  between 
curves  moves  to  smaller  L.  This  is  simply  explained  by  the 
window  for  the  phase  variation.  For  large  L,  the  value  of  kp 
at  which  the  window  decreases  is  on  the  order  of  k^,  the  location 
of  the  peak  of  kpEnn(kp).  As  L  is  decreased,  the  window  affects 
the  function  at  larger  values  of  kp;  therefore,  the  effect  occurs 
at  smaller  values  of  p.  Once  the  window  decrease  occurs  for 
kp  >  k5,  the  effect  is  lost.  Figure  11  further  illustrates  this 
by  showing  the  structure  function  at  small  p  as  a  function  of  L. 

It  is  significant  that  the  change  occurs  at  physically  meaningful 
values  of  L. 

The  L  dependence  of  the  shape  of  the  logarithmic  amplitude 
structure  functions  is  stronger  as  the  logarithmic  amplitude 
window  has  a  stronger  L  dependence.  As  L  is  decreased,  increas¬ 
ingly  more  of  the  lower  wave  number  (longer  length)  portion  of 
the  spectrum  is  eliminated  or  "cut  off"  by  the  window.  This  is 
apparent  as  long  wavelencrth  phenomena  cannot  affect  short  distance 
measurements,  but  will  appear  as  nonstationari ties  in  the  mean 
and  be  eliminated  in  a  structure  function  measurement  (see 
appendix  A) .  Thus,  as  L  is  increased,  the  lowest  wave  number  at 
which  the  integrand  provides  "energy"  is  increased.  Hence,  the 
separation  distance,  p,  at  which  the  structure  function  flattens 
is  decreased  (see  figures  7  and  8) .  As  the  "cutoff"  of  the  window 
is  given  approximately  as  ~\fTT/K,  the  "cutoff"  value  of  p  for 
IA5  <  p  <  1/k^  iS\fc/h.  For  L  >  10(c/k|,  the  window  no  longer 
changes  the  integrand.  Similarly,  for  L  <  .  lic/kj*,  the  spectrum 
itself  "cuts  off,"  so  the  window  no  longer  has  an  effect  i  ther 
than  scaling. 

Figure  10  shows  the  dependence  of  the  variance  on  L.  For 
L  >  K/kJ,  L  is  seen  to  linearly  scale  the  variance.  For  K/k^  <  L 
<  K/k?,  as  Tatarski2  predicts,  the  variance  varies  as  L11'  .  For 
L  <  K/k|,  the  variance  changes  as  L  .  Once  again,  a  rise  cor¬ 
responding  to  the  viscous-convective  range  is  evident.  For 
K/kg  <  L  <  K/k 3,  the  variance  is  above  the  and  L3  asymptotes. 

This  is  significant,  as  in  physically  .iteaningful  ranges  of  L, 
from  1  centimeter  to  100  meters,  the  variance  of  the  salinity 
variations,  for  example,  is  above  the  asymptotes  by  a  factor  of  10. 

CONCLUSIONS 

©  In  sea  water,  the  existence  of  the  viscous-convective 
subrange  in  the  temperature  and  salinity  (hence  refractive  index) 
spectra  causes  a  significant  increase  in  the  fluctuations  of  the 
phase  and  logarithmic  amplitude  of  the  optical  wave.  The  increased 
optical  fluctuations  caused  by  the  salinity  microstructure  occur 
for  optical  path  lengths  of  about  1  centimeter  to  1000  meters; 
thus  experimental  confirmation  is  feasible. 
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©  The  use  of  structure  functions  instead  of  correlation 
functions  to  characterize  the  spatial  fluctuations  is  considered 
a  more  meaningful  approach  because  the  physical  phenomena  being 
described  are  at  best  only  quasi-stationary ,  with  stationary  first 
increments  assumed.  The  structure-,  function  approach  directly 
provides  the  functional  dependences,  particularly  at  small  distances, 
because  these  dependences  are  not  hidden  by  the  variance  as  they 
would  be  in  the  case  of  correlation  functions. 

©  The  smooth,  mathematical  representation  developed  for 
the  refractive  index  spectra  satisfies  theoretical  and  experimental 
data  where  available  and  is  suitable  foi  computer  solution  of  the 
integral  equations.  Further,  this  spectra  representation  is 
sufficiently  general  so  that  it  can  be  used  for  various  media  by 
appropriate  selection  of  the  spectral  cutoff  and  breakpoint  wave 
numbers. 
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F'c;ure  2 

Three-Dimensional  Spectrum  of  Refractive  Index 
Variations,  Enn(k),  as  d  Function  of  Wave  Number,  k 


3-C 
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Figure  3 

Modified  Three-Dimensional  Spectrum  of  Refractive 
Index  Variations ,  En(k)f  as  a  Function  of  Wa -e  Number,  k 
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Figure  4 

Structure  Function  of  Refractive  Index  Variations, 
Dni^r)»  as  a  Function  of  Distance,  r 


figure  5 

Computed  Phase  Structure  Function  for 
Temperature,  as  a  Function  of  Separation 

Distance,  p,  and  Distance  Traversed,  L 
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Figure  9 

Length  Normalised  Optical  Phase  Variance 
as  a  Function  of  Distance  Traversed, 
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Figure  11 

Optical  Phase  Structure  Function,  D^#  at  Separation 
Distance  p  =  0.2  x  10  **  as  a  Function  of  Distance 

Traversed.  I, 


APPENDIX  A 


STATISTICAL  DESCRIPTIONS  OF  RANDOM 
FUNCTIONS  AND  RANDOM  FIELDS 
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LIST  OF  SYMBOLS 


( • ) *  Indicates  complex 
conjugate 

(• )  Indicates  average 

x ,  t  Time  variables 

a)  Radian  frequency 

ip  ( - )  Fourier  transform 
of  time  function 

6  ( • )  Dirac  delta  function 

i  ~\f~l 

(f(t)}  Ensemble  of 

functions,  f (t) 

m  ( • )  Mean 

B  ( • )  Correlation  function 


C(-) 
D  { •  ) 


Bf  (•) 
Bff(* 


Covariance  function 

Structure  function 

Vector  distance 

Vector  wave  number 

Modified  three- 
dimensional  spectrum 

)  Three-dimensional 
spectrum 

One-dimensional 

spectrum 

Two-dimensional 

spectrum 


ANALYSIS  OF  A  RANDOM  FUNCTION 

A  set  or  ensemble  of  random  functions,  {fp(t)},  can  be 
described  by  all  possible  multidimensional  probability  distribu¬ 
tions.  In  practice,  only  the  mean  and  the  autocorrelatio  func¬ 
tion  (the  first  joint  moment)  are  employed.  If  the  Gauss  .an 
distribution  applies,  the  functions  ere  completely  descr  oed. 

A  priori  knowledge  or  assumptions  on  the  properties  of  tue  random 
function  provide  different  methods  for  obtaining  the  descriptive 
functions. 

The  mean  and  the  autocorrelation  function  are  given  by  the 
expressions  in  (A-l) .  The  averages  here  are 


mf(t)  =  ifp (t) 
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A-l 


(A-X) 


Bf(t,T)  =  {f  (t)f*(t  +  T)} 

where 

■f*  V> 

f  C t )  =  the  p  realization  of  the  random  function 

^  at  argument  t 

m^(t)  =  the  ensemble  average  at  argument  t 

B£(t,x)  =  the  autocorrelation  function  at  argument 

t  and  log  x. 


If  mj(t)  and  Bf(t,x)  are  independent  of  t,  the  ensemble  is  said 
to  be  weakly  stationary.  If  all  the  probability  distributions 
are  independent  of  t,  the  ensemble  is  stationary.  If  the  ensemble 
is  guasi-erodic,  then  the  mean  and  autocorrelation  functions  may 
be  obtained  by  the  expressions  in  (A-2)  employing  averaging 
across  any  member  of  the  ensemble. 


m*;  =  lim  — 
r  t+co  T 


Bf(x)  =  lim  |  j'  f(t)f*(t  +  x) dt 
-T/2 


(A-2) 


If  all  the  probability  distributions  can  be  obtained  by  such  an 
average,  the  ensemble  is  ergodic.  If  the  properties  can  be 
derived  from  any  short-term  average  (removing  the  limits  of  the 
expressions  in  (A-2)),  the  ensemble  is  referred  to  as  weakly  or 
strongly  self-stationary,  where  the  weak  connotates  the  same 
meaning  as  in  the  case  of  stationarity . 

ANALYSIS  TECHNIQUES  FOR  WEAKLY  SELF-STATIONARY  FUNCTIONS 

Where  quasi-ergodicity  and  weak  self-stationarity  can  be 
assumed,  additional  descriptions  of  the  function,  based  on  the 
mean  and  autocorrelation  function  are  available.  Expression 
(A-3)  gives  the  variance;  expression  (A-4)  gives  the  covariance 
function.  In  the  case  of  a  mean  free  function,  the  autocorrelation 
function  and  the  covariance  function  are  equivalent. 
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A-2 


c£  (0)  =  Bf (0)  -  m\ 


(A- 3) 


Cf  (T) 


Bf(t) 


-  m.  = 


£  /  (f(t) 


mf ) (f (t  +t) 


m^) dt . 


(A- 4) 


A  stationary  function  may  be  represented  in  terms  of  a 
Fourier-Stieltjes  integral  as  shown  in  expression  (A-5) . 


f  (t) 


CO 


—  00 


(A-5) 


The  autocorrelation  function  can  then  be  written  as  expression 
(A- 6) 


iio,  t, 

Bf(tl“  t2)  “  f(t1)f*(t2)  =/ e  1  xdt^ 


-itu^t, 

e  ~  (oj2  ) 


a>2^'2^ 


dij/f  (o)1)di|;|  (w2) 


(A-6) 


For  stationary  ensembles,  the  ensemble  average. 


d4>  (u1)dti/*  ((d2)  , 


can  be  written  in  terms  of  a  power  spectra  Wf(w^),  as  expression 
(A- 7)  . 


d^ (w^d*^ (w2)  -  Wf(u1,6(u1  -  a)2)dd)1da)2.  (A-7) 
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In  the  case  of  weak  self-stationarity,  BCt-^  -  t2>  must  depend  only 
on  the  difference  x  =  t^  -  t^,  so  that  the  relationship  (A-7)  can 
be  employed.  Thus,  the  Fourier  transform  pair  relationship  between 
the  autocorrelation  function  ana  the  power  spectra,  the  expressions 
in  (A-8)  can  be  written.  As  both  functions  are  even,  the  expres¬ 
sions  in  (A-9)  follow. 


B^(x)  =  ^e^WTWg  (ui)dw 


o 

wf (w)  =  h  J~ 


elU>TBf  (x)dx 


(A-8) 


Bf(x) 


Wf  (w) 


-/■ 


cos(asx)W£(aj)daj 


4 


cos  (iox)Bf  (x)dx. 


(A-9) 


ANALYSIS  TECHNIQUES  FOR  A  CLASS  OF  NONS TAT T.ONARY  RANDOM  FUNCTIONS 

Many  random  functions  descriptive  of  physical  phenomena 
cannot  be  assumed  to  exhibit  even  weak  stationarity .  For  example, 
the  temperature  and  salinity  of  the  ocean  are  characterized  by 
slowly  changing  means  and  variances.  A  technique  has  been  devel¬ 
oped  whereby  a  meaningful  analysis  can  be  obtained  for  the  class 
of  random  functions  known  to  have  stationary  increments.  This 
will  be  defined. 

In  the  general  case  of  a  nonstationary  random  function,  the 
autocorrelation  is  defined  (reference  (a) )  by  expression  (A-10) . 


BA(tl,t2)  =  f(t1)f*(t2)  = 


i 


“2t2) 


d\|/f  (u1)d^*(ui2) 
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A- 4 


00 


SA(fl'f2)dfldf2'  (A— 10) 


It  is  desirable  here  to  make  the  variable  changes  defined 
by  the  expression  in  (A-ll)  which  provide  the  definitions  of  the 
expressions  in  (A-12) . 


T 


+  t. 


\ 


t 


1 


=  t 


-  T/2 


f 


t,  =  t  +  t/2 


g  = 


f  1  +  f  2 


(A-ri) 


f2  =  g  +  f/2  ) 


BB(x,t)  =  BA(tx,t2) 

sB(f,g)  =  sA(frf2)  j 


(A-12) 


The  expressions  in  (A-10)  and  (A-12)  provide  the  expressions 
in  (A-13) . 


BB(T 


i2n  (fx  -  gtj 


Sfi(f ,g)dfdg 


•# g)  =  e 


SB(f,g)  =1  le-i2,r(fT  - 


Bg (t,T)dxdt. 


(A-13) 
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If  f  (t)  were  stationary,  the  expressions  in  (A-14)  must 
hold  to  satisfy  expression  (A-8) . 


Bf(T}  =  Bg  (t  ,-t)  S  (t) 
Wf(f)  =  Sg(f,g)6(g). 


(A-14) 


The  Dirac  delta  functions  of  the  expressions  in  (A-14)  show 
the  lack  of  dependence  on  ti  and  t2.  If  this  lack  of  dependence 
occurs  only  over  a  range  x  <  Xg,  the  expressions  in  (A-12)  can 
be  approximated  (reference  (b) )  as  shown  by  expressions  in  (A-15) . 


B3(r/t)  =  Bl(T)B2(t) ;  x  <  xQ 
SB(f,g)  =  S1(f)S2(g);  f  <  l/xo 


(A-15) 


The  applicability  of  this  approximation  can  be  demonstrated. 
This  leads  to  a  practical  method  for  obtaining  stable  information 
about  the  random  function  in  question. 

Consider  a  mean-free ,  random  function  with  a  slowly  varying 
multiplicative  scale  factor.  The  variance,  for  example,  is 
varying  with  this  factor.  The  stationary  variations  can  be 
separated  if  the  scale  length  of  the  variations  of  the  scale 
factor  is  x  >  xQ,  where  xQ  >  x^,  the  scale  length  which  includes 
the  stationary  variations.  A  correlation  of  the  form  given  in 
expression  (A-10)  performed  at  some  t  =  ta  will  have  the  same 
form  as  one  performed  at  t  =  t^  except  for  a  scale  factor.  All 
such  correlations  will  be  0  for  Xi  <  x  <  x  .  Af  x  >  xQ,  the 
approximation  no  longer  holds,  ana  a  deviation  from  zero  occurs. 
Thus>  the  correlation  function  of  the  stationary  variations  can 
be  obtained,  and  the  scale  xength  of  the  nonstationarity  can  be 
determined. 

Difficulty  occurs  if  the  random  function  contains  a  non^ 
stationary  mean.  Correlations  of  the  form  of  expression  (A-10) 
cannot  be  usefully  employed.  The  average  of  expression  (A-l8) 
provides  a  method  of  analysis. 
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[*<*!>  -  f(t; 


12 


BB(0,t  -  T/2)  +  BB(0,t  +  T/2) 


-2BB(T,t). 


(A-16) 


An  important  feature  of  the  average  of  expression  (A-16)  is 
that  any  variations  in  the  mean  of  scale  length  greater  than 
t  =  t2  -  tj_  will  not  affect  the  result.  Thus,  expression  (A-16) 
can  be  treated  as  if  the  random  function,  f (t) ,  were  mean  free 
and  expression  (A-17)  can  be  written  for  r  less  than  the  scale 
length  of  the  variations  of  the  mean. 


jf  (tx)  -  f  (t2)j2  =  B^O)  B2(t  -  t/2)  +  S2(t  +  t/2) 


-2B1(T)B2(t) 


(A-17) 


Expression  (A-18)  can  be  written 

E2(t  -  t/2)  =  B2(t)  =  B2(t  +  t/2)  ;  t  <  TQ  (A-18) 


if  the  scale  length  of  the  variations  of  the  scale  factor  exceeds 
t.  Finally,  expression  (A-19)  can  be  written  if  the  scale  length 
of  the  variations  of  the  mean  and  the  scale  factor  (variance) 
exceed  vQ. 


f  (tx) 


-  f  (t_2)|  =  2B2  (t)  jj 


2B2(t)  Bx(0)  -  Bx(t) 


]■ 


T  <  T. 


(A-19) 


If  the  average  of  expression  (A-19)  is  evaluated  at  different 
t's,  the  result  will  remain  t  invariant.  If  the  scale  length  of 
the  stationary  variations  is  t^  <  tq,  the  result  will  be  zero  at 
t  =  0,  a  constant  for  t^  <  t  <  tq,  and  will  deviate  from  the 
constant  for  t  >  t0.  Generally,  a  random  function  which  can  be 
analyzed  in  this  manner  is  said  to  have  stationary  increments. 
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The  average  introduced  by  expression  (A-16)  is  referred  to 
as  the  structure  function  (reference  (c) ) .  The  expressions  in 
(A-20)  define  this  function  using  the  definitions  of  the  expres 
sions  in  (A-ll) . 


DA(tit2)  = 


'B 


(T,t)  =  jf(tr)  -  f(t2) 


12 


=  Df (t,T) . 


(A-20) 


SUMMARY  OF  MOST  OFTEN  EMPLOYED  ANALYSIS  TECHNIQUES 

The  most  often  employed  analysis  techniques  are  summarized 
in  this  section.  It  vi3 1  be  assumed  that  the  random  function  is 
at  least  weakly  self-stationary  except  for  nonstationariti.es  in 
mean  and  a  scale  factor  (variance) .  Quasi-ergodicity  is  then 
also  assumed. 

Mean  and  Variance 

The  mean  is  calculated  as  shown  in  expression  (A-21) . 


(A-21) 


T  is  chosen  based  on  the  scale  length  of  the  period  of  the 
lowest  "sationary*'  variation. 


The  variance  is  calculated  as  shown  in  expression  (A-22) . 


C*(0) 

Jm  ' 


(t)dt  -  m|. 


0 


(A-22) 


T  must  be  chosen  as  expressed  above. 


A-8 
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Autocorrelation  and  Covariance  Functions 

The  autocorrelation  function  is  calculated  as  shown  in 
expression  (A-23) . 


Bf(r)  = 


T 

=  |Jf(t)f*(t 


+  x)dt. 


(A-23) 


Again  T  must  be  carefully  chosen.  The  covariance  function  is 
calculated  as  shown  by  the  expression  (A-24) . 


T 

cf(T)  =  (f  (t)  -  m) 


(f(t  +  t)  -  m)*dt 


=  Bf (x)  -  n  . 


(A-24) 


Thus,  for  mean- free  random  functions,  the  correlation  and  covariance 
are  equivalent.  In  the  case  of  nonstationarities  in  mean  and 
scale  factor  (variance) ,  particularly  if  the  mean  is  large  as 
compared  to  the  variance,  expressions  (A-22) ,  (A-23) ,  and  (A-24) 
may  not  provide  useful  results. 

Structure  Function 

Expression  (A-25)  defines  the  structure  function. 


5 

Di(T)  -  ij 


[f(t)  -  f(t  +  T)J  Ct. 


(A-25) 


T  is  determined  as  expressed  above.  The  structure  function 
to  be  most  often  employed  is  defined  by  the  expressions  in  (A-26) , 
where  t  is  the  scale  length  described  above. 


D,(t)  =  D*{x);  t  <. 


(t)  —  ;  x  > 


(A-26) 
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The  definitions  of  the  expressions  in  (A-26)  provide  a 
methodology  of  defining  a  mean-free  autocorrelation  function  in 
teruir  of  the  structure  function  which  can  be  universally  applied 
in  later  analysis.  This  is  shown  in  the  expressions  in  (A-27) . 


Pf  (t)  =  2[Bf  (0)  -  Bf  (t ) ) 

Df(c)  =  2Bf (0)  =  2Cf(0)  (A-27) 

Bf(T)  =  1/2 (Df (a)  -  Df (T)l 


Spectral  Functions 

The  expressions  in  (A-28)  were  given  above  as  expressions 
(A-9)  . 


C 

=  2/c 


Bf(r)  =  2 j  cos(ex)lf^(c3)d£3 


Kf(u)  =  ||^COS(uT)Bf  (x)dT. 
0 


(A-28) 


The  expressions  in  (A-29)  relating  the  structure  function 
to  the  spectral  function,  follow  from  expression  (A-28)  and  the 
definitions  of  expressions  (A-27) . 


Df(T)  =  4 J~ a  - 


cos  (et)  )W^(a)da 


D*(a) 

«£<“>  -  -T- 


-  £/« '<* 


(A-29) 


©(y)  -  -=rr  I  cps(«T)D^(T)dT, 
0 


Report  6-196 


A— 10 


ANALYSIS  OF  RANDOM  FIELDS 

A  random  function  of  three  variables  is  a  random  field. 
Only  scalar  fields  are  considered  here.  A  random  field  is 
'.on  side  red  homogeneous  if  its  mean  value  is  constant  and  auto¬ 
correlation  does  not  change  with  position.  The  expressions  in 
(A-30)  show  these  results. 


f  (r)  =  constant 


Bf(rl,r2)  =  B 


f  «i  +  V 


+  ro> 


(A-30) 


Setting  r  =  — r2  in  expression  (A-30)  provides  the  expressions 
in  (A— 31) .  ° 


B£(r^,r2)  0)  =  ir-,  —  r^)  —  B^(r).  (A— 31) 


'f  %tl 


- 2 ' 


f  VJ-1 


Expression  (A-31)  shows  that  homogeneity  corresponds  to  weak 
stationarity  in  random  functions. 

If  Bf(r)  depends  only  on  |rl  =  r,  the  random  field  is  iso¬ 
tropic  .  Thus,  isGtropy  corresponds  to  quasi-ergodicity . 

A  random  field  nay  be  represented  by  a  Fourier-Stieltjes 
integral  as  shown  by  expression  (A-32) . 


The  autocorrelation  function  can  then  be  written  as  shown 
by  the  expressions  in  (A-33) . 


f(r) 


■Iff 


exp (ik*r) d$f (k) 


(A-32 


Bf(r)  =  f  (r^,)  f*  (r2) 


-  3c2*r2ld^f  (kx)d^|(k2) 


A-iX 
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(A— 33) 


exp  (ik>r)  (k)  dk 


where 


d$f  (k^d^Xkj)  =  6{k±  ~  k2)Eff  (k1)dk1dk2. 


Considering  the  even  characteristic  of  the  functions,  the 
expressions  in  (A-34)  show  the  Fourier  transform  pair  relationship 
between  the  autocorrelation  and  the  three-dimensional  spectrum. 


Bf(r)  = 


OQ 

= JJJcos(k‘ 


r)Eff (k)dk 


'ff  ’ 


J 


E^clk)  =  1^-1  !  j | cos (k • r ) Bf (r ) dr . 


(A-34) 


ISOTROPIC  HOMOGENEOUS  RANDOM  FIELDS 

If  the  random  field  is  isotropic  and  homogeneous  on  a  plane 
x  =  const,  the  autocorrelation  can  be  expressed  in  terms  of  a 
two-dimensional  spectrum  which  can  be  related  to  the  three- 
dimensional  spectrum.  This  is  developed  in  the  expressions  in 
(A-35) . 


B, 


oo  as  2l 

<*.P>  -J" j f, 


d9kpdkpdkx  cos(kxx  +  kpp  cos6)Eff (kp,kx) 


-«*  0  0 

ee  00  2a 


■Iff 

-00  0  0 


d3k  dk  dk - [cos (k  x)  cos (k  pcosS) 

H  K  A  H 

-  sin(kxx)  sin(pkp  cos8) ]Eff (kp,kx) 
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(As  Eff  is  even) 
eo  oo  2n 


d6kpdkpdkxcos (kxx)  cos (pkp  coed ) Ef f (kp ,kx) 


-«*  o  0 


(As  2irJ 


2n 

o(kpp)  =  j cos(kj 


P  cos0)d8 


2.  Hops  (kxx)  Jp  {KpP )  Ef  f  ;<kpkx)  kpakpdkJ 

-W  0 


to  oo 

#  2rj" “'pVo'V’  JdkxE ff<V 

0  -oo 

a> 

=  dkpkpJ0(kpp)Ef2(kp,x). 


V  cos-(kxx) 


(A-35) 


The  final  expression  shows  the  two  dimensional  spectrum,  E,~ 
(kp,x).  The  Fourier  transform  relationship  to  the  three-dimen¬ 
sional  spectrum  is  shown  in  expression  (A-36) . 


Ef 2 (fcpX) 


a 

f 


=  /  cos  (k  x)E,^(k  ,k  }dk 
/  x  f-f  xf  p  3 


Eff (kx,kp) 


a 


Ef2(kp,x) 


cos  (k^xjdx. 


(A-36) 
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Further  insight  into  the  two-dimensional  spectrum  is 
shown  by  the  expressions  in  (A-37) . 


f (x,y/Z)  =  /  jexp[i(kyy  +  kzz) ]d  f2(ky/kz/x) 


6  (k2  -  k:2)6  (k3  -  k^)Ef2  {ky/kz,x)dkydkzdk^dk^ 


=  d*f2(ky»kz,x)di|»*2(k;^k^,x) . 


(A-37) 


The  first  of  these  shows  the  fandom  field  on  the  plome  x  =  const 
expanded  in  a  Fourier-Stielt jes  integral.  The  second  expression 
shows  the  two-dimensional  spectrum  in  terms  of  that  expansion. 

If  the  random  field  is  homogeneous,  the  correlation  function 
and  the  two-dimensional  spectrum  are  independent  of  x  as  shown 
in  the  expressions  in  (A-38) . 


Bf ( P )  =  2*l  dkpkpJo(kpp)Ef2(kp) 
0 


00 

r 


Ef2(kp)  _  I  Eff (kx,kp)dkx* 


-po 


(A-38) 


If  the  fandom  field  is  homogeneous  and  isotropic,  the  auto¬ 
correlation  can  be  expressed  in  terms  of  a  one-dimensional  power' 
spectrum  which  can  be  expfessed  in  terms  of  the  three-dimensional  ■■ 
spectrum.  This  is  developed  in  the  expressions  in  (A-39) , 
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00  ir  tt/2 


Bff(kr) 


cos(krr  cos0 ) Bf (r) (r  sinSdO) (rd0)dr 


_oo  o  — ir/2 


2  /  cos(krr  cps0)Bf  (r)r2  sin0d0dr 


2(2u) 


00  7T 

2  J"dr  ^  fB£(r)M  ~sin  (kr 


r  cos0)  d6 


oo 

-^rjfdr 
(2tt  )  2J  kz 


j—  Bf  (r)  sin  (k  r j 
r  r 


1  1  dll 

3? TT  HF  dr  cos  <krr )  Bf  (r) 


=  -1-  2 

2lTkr  2kr  Efl(kr} 


00 


rB|  (r)  sin^rjdkj 


OO 

Bf(r)  =  TrJ"  krEff  (kr}  sin(kf 


r)dkr. 


(A-39) 


The  expression?  in  (A— 39)  show  that  one— dimensional  spectrum 
exhibits  the  expected  Fourier  transform  relationships  of  the 
expressions  in  (A-40) . 
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00 


Efi'V 


dr  cos (krr)Bf (r) 


—00 


CO 


Bf  (r) 


cos  (krr)Efl (kr) . 


—00 


(A-4C) 


In  addition,  expression  (A-41)  may  be  written 


Efi  {kr)  =  /  dk  (2uk)  Eff  (k)  +  /  dk  (2irk)Eff  (k) 


-k 


=  |  dk(47ik)Eff  (k)  . 


(A-41) 


All  the  results  of  previous  sections  defining  relationships 
between  the  autocorrelation  and  power  spectral  functions  apply 
to  the  autocorrelation,  Bf(r),  and  the  one-dimensional  spectrum 
Efl(r)  . 

LOCALLY  HOMOGENEOUS  AND  ISOTROPIC  RANDOM  FIELDS 

The  methodology  of  previous  sections  on  local  stationarity 
is  applied  to  locally  homogeneous  random  fields.  The  three- 
dimensional  structure  function  is  defined  by  expression  (A-42) . 


-Df(Ej_,r2)  =  [f(rx)  -  f(r2))^. 


(A-42) 


Report  6-196 


A-16 


The  field  is  locally  homogeneous  if  shifts  of  and  t2 
within  a  region  do  not  change  the  structure  function  as  shown 
in  expression  (A-43) . 


Df(rl  +  V  r2 


+  r  )  =  (r.  ,r0)  =  D^fr.  -  r0) . 


(A-43) 


In  case  of  local  homogeneity  and  isotropy  on  a  plane  x  = 
const,  expression  (A-44)  results. 


00 

Df(p)  »  J  a  -  Jo(kpP)]E£2(kp,0)kpdkp. 


(A-44) 


In  the  case  of  local  homogeneity  and  isotropy,  expression 
(A-45)  results. 


Df(r)  =  2  (1  -  coskrr)Efl (kr)dkr 


00 

f  r  sink  tl  2 

■  4*J  L1  -  E£f(V*rdlV 


(A-45) 


The  one-,  two-,  and  three-dimensional  spectra  are  related  as 
above  except  that  no  Dirac  delta  exists  for  the  mean.  The 
structure  and  correlation  functions  are  considered  modified  as 
expressed  above  for  large  p  and  r. 

2  v 

The  term  47rkrEff(k  )  or,  in  general/  47Tk”£ff(k)  seen  in 

the  last  of  expression  TA-45) ,  is  employed  often.  For  simplicity. 
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APPENDIX  B 

COMPUTER  PROGRAM  FOR  C/ujOULATTNG 
THE  OPTICAL  PROPAGATION  PROPERTIES 


This  program  is  utilized  for  all  four  of  the  required 
calculations.  A  particular  calculation  is  obtained  by  changing 
lines  17  and  19  in  the  program.  In  line  17:  ICODE  =  -1, 
logarithmic  amplitude  variations  are  calculated;  if  ICODE  =  +1, 
phase  variations.  Line  19  sets  the  kg  factor;  thus  XKS  =  18  for 
temperature;  XKS  =  180  for  salinity. 

The  program  provides  both  printed  results  and  a  control  tape 
for  a  plotter  to  provide  plotted  results. 
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CAPwOPT«T3d.lO.TPl.P7. 

CHARGE.CARM,184?OOAOI  .PR.R#  842’* 

ATT  ACH  (  SCOPS ♦ CASGSCOPS . 10=24000 1 202) 

PFOUFST .TAPE4R.HI »S. (CAD 343/P INGIN) 

FTN(T) 

HEDUCF. 

LOAO (SCOPS) 

LGO. 

UNLOAD (TAP£4 A) 

OOOOOOOOOOOOOOOOOftOOOO 

CoXr,MOEtr«5;»™PE6=OU’-’,|T-T‘?E‘n’'NPUT-T*PFS'"'Pl''' 

OIMCNSTON  W(9S) .XRHO(A5) .EE{4) .00(4) .TITLF(A) 

CALL  CAMOAV (35)  L 

P£A0(5« 10)  (TITLE(I) *1=1 .4) 

10  FOPHAT(AAlD) 

XMlN=-5. 

XMA x=4. 

Y»*IM=2, 

VMAX=30. 

nx=l. 

DY=1. 

CALL  SET*Iv(24,0, 100.24) 

PI=3.141S02 

TCODF=i 

8M*X=1000. 

Xf.5=130. 

I  J*=l 
*  >“fll  =  1.8 
WPTTE(S.IOI) 

101  FO°MAT(lhl) 

WRITE (6. 102) 

102  P0=MAT(1m  .*  L  RHO  OPHISHO.) 

WRITER. 103) 

103  FORMAT ( lh  ) 

XL=1 .E-2 

00  5  15=1.16 
JJ=0 

FXO0=! .E-5 
00  5  16=1.11 
no  7  17=1 .s 

PH0=2.*I7*ExP0 

XPHO(JJ)=RhD 

oon>=.ooi4 

suM=n. 

A=0. 

B=A*00(IjK) 

47  AA=.S«(A.B) 

48=8-4 

CC=. 4301449*89 

Y=. 05061427* (FUNC(AA«CC) *FI)NC  ( AA-CC) ) 

CC=. 3943332*98 

Y=Y* .III  1 90S* (FUNC ( *A«CC) *FiJNC (A A-CC)  ) 

CC=. 2627662*89 

Y=Y*,1S69533*(F(1NCCAA.CC)  *F|»nC(AA-CC)  ) 

CC=.09l 7 1732*88 

Y=R3* (Y+, J 4)3419* (FUNC(44*CC>*FiJNC{AA-CC)) ) 


842.KY8RANIEC 
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SUH=SUM*Y 

IFC9-«M*x)  46.46,49 
46  DOCIJK)=DOCIJK)*EE  CIJK) 

A=R 

9=A»D0(tJK) 

60  TO  47 

4«  W(JJ)=4.*PI*SIW 

VPITE(6«4)  XL«»HO,WCJJ> 

4  FORMATClH  *£9.l .2X.E9.1 .2X.E12.5) 

7  CONTINUE 

EXP0=EXP0»10. 

6  CONTINUE 
WPITEI6.104) 

104  FORMAT (1H  ) 

XBHO  C 1 ) = AL0G1 0 C  XRHO  Cl)) 

MCI) =AL06 Ifl(WCl)) 
no  11  111=1,54 
112=111*1 

XBHO  C 1 12) =4L0G1 OCXRHO C 1 12) ) 

WCIC2)=AL0010  CWCI12) ) 

IX'i=lXVCXPHOCIl  1)  ) 

IX2=IXVCXRH0CI121) 

IY1=IYVCNCI11)) 

IY2=lYVCWCI12) ) 

CALL  LIMEVCIX1.IY1.IX2.IY2) 

11  CONTINUE 
XL*XL*10. 

5  CONTINUE 
CALL  PLTNOCO) 

STOP 

END 

FUNCTION  FUNCCXKP) 

COMMON  ICOOE-XL.XK5.RHO 

XK1=.001 

XK3=.2« 

XK=.n?161F6 
PI=3. 141592 
Z=176. 

Ol  =— 2 • *  C  XKP/XKS) **2 
IFCZ-AHSCOl))  32.31.31 
32  ENN=0. 

60  TO  34 

31  FNN=C  CXK3**C2«/1.)/XKP*»CI 1 ,/3.) l*Cl./f  Cl •*CX<l/XKP)**2),*Cll./6.) 

1) -l ./C  Cl.*CXK3/XKP)**2)«*ni./6.) ))*Cl./CCXKP**2»XK3**2)**C3./2.)) 

2) )*EXPC01) 

34  X=XKP**2»XL/XK 

CALL  CALSINCX.iMIWjS) 

TFCICOOE)  ’1. 21 .22 
21  E2=?I*CXK#*2)*XL*amINUS*E*p< 

60  TO  23 

2?  APLUS=2. -AMINtlS 

E2=PI *  C  XK*  »2) *XL*APLUS*ENN 
23  ARG=xkp*rho 

CALL  3ESCAPG.SUK) 

Fl|NC=  C 1  .-SUM)  *F2*X<» 

IFCABG.LT..1)  FlJNC=SUM*F2*XKP 

BETU»N 

END 

SIW BOUT INF  CALSINCX.AMINUS) 

XCONl*l.F,-fio 
91=1 .—X 
IFCRl)  6,6,9 
9  K=  1 

AMlNIJSsO 
4  FACT=1 

T£MB=C  C-l .) **CK*1 ) )*CX**C2.*K) ) 
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M=*MS'(  ltHv> 

XCON=I.E-50 

*2=xcon-ai 

ir!M2)  11,7.7 

11  JJ=2*K*1 

f>0  9  11=1 .  jj 
*J=IT 

FACT=F*CT*Aj 
ft  CONTINUE 

TEMS>=TFNP/F*CT 
TFST=XCOUl-*BS(TEMP) 
IFITEST)  12.7.7 

12  KrK»l 


* 

7 


GO  TO  4 


*Mi»ius=n.-tsiMcxr.'xii 

RETURN 
EOT 

SUBROUTINE  BESfX.SUM) 
CONUl.E-50 
CON=l.f-S0 
PI=3. 1*1592 
IFCX-IO.SI  1.1,2 
1  K=l 
SUM=! . 

IFCX.LT.. 11  sofci). 

3  F*CT=1. 

XWMr  (  f-1 .  )  *«}  ,  j  ;.25,Xf  ,2,  , 

IFfX.LT^.lJ  XNUM=-xwt 

*I=*«JSCXNI»*l 

*2=C0N-*1 

TF142I  *.ft,8 

*  no  6  I f =1 ,K 

*J=II 

F*C7=F*CT**J 

6  CONTINUE 

tf »s»w/  (fact**?.  > 
TEST=C0NI-*8SCTE-HPJ 
TFITESTl  7.8.8 

7  K~K*t 
SU“=S'JH.TFyp 
9=F*CT**2. 

IFC1.F65-OJ  8.8.50 

59  IFC1.E37— XfftWj  8.51,51 
51  GO  TO  3 
2  D0=l. 

K=1 

15  **=*«1C 
**.*1=1 

no  9  L=I.*f*.2 
NU*<l=Ni»«I*L**2 
9  CONTINUE 

*OTMl  =  (-l.«*K)*Nyp(j 

*1=*«S(XNUMI1 
*R=*1-!.F85 
IFCA21  10.11,11 
19  F*CT=1. 

II=?*K 

no  12  11=1, JJ 
*J=II 

F*CT=F*CT*AJ 
12  COHTINUT 

TEST=C^II-*55CTF**»l 
IFCTEST1  *3.11.11 
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13  lCs<»I 
GO  TO  15 
11  0O*-l./f8.*XI 

K±Z 

16  NU»*2= 1 
K3=**K-2 

OO  17  L-1*K3*2 

NW2*W2«t»*2 

17  COMTIMjE 

xmr*z=  (-1  .**k»  •mm? 

*i=ABscxmm2» 

A2=A1-1.E65 
in*?)  18,19.19 

18  FACT=1 

00  20  IIxl.jj 
AJ*II 

FlCTsFUCfftJ 
?0  CONTIHUE 

TEm»*mr«/iF*cT*  i8.*xi  **  i2.*«-i  1 1 

QO=90*TE*P 
T£ST=C0N1-A8S I TEMPI 
IFfTESTI  21,19.19 
21  KsX.l 
GO  TO  16 

i«  stm=ii?./fPi*xii«*„5>*iPO*coscx-»»i/*.i-ao*siMfx-*>i/*.j) 

8  BFTl|R«l 

E'T) 

0000000000000000003000 

LOGCOPH1 tRMOl }  VS-  LOGISHO)  F0«  S*t_.  AMO  L  =  10**-2  TO  lfc**13  BY  FOftCRS  OF  10 

0000000000000000000000 
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